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Data from the first physics run at the Relativistic Heavy-Ion Collider at Brookhaven National 
Laboratory, Au+Au collisions at ^/snn = 130 GeV, have been analyzed by the STAR Collaboration 
using three-pion correlations with charged pions to study whether pions are emitted independently 
at freezeout. We have made a high-statistics measurement of the three-pion correlation function 
and calculated the normalized three-particle correlator to obtain a quantitative measurement of the 
degree of chaoticity of the pion source. It is found that the degree of chaoticity seems to increase 
with increasing particle multiplicity. 

PACS numbers: 25.75.Gz, 25.75.Ld 



Two-pion Hanbury Brown and Twiss (HBT) interfer- 
ometry in principle provides a means of extracting the 
space-time evolution of the pion source produced at kine- 
matic freeze-out in relativistic heavy-ion collisions Q, Q . 
An underlying assumption of this method is that pions 
are produced from a completely chaotic source, i.e. a 
source in which the hadronized pions are created with 
random quantum particle production phases. In applica- 
tions of two-pion HBT the validity of this assumption is 
usually tested by extracting the "A-parametcr" which in a 
simple picture is unity for a fully chaotic source and zero 
for a fully coherent source 0,0- However, this parameter 



also depends on many other factors, such as contamina- 
tion from other particles in the pion sample, unresolvable 
contributions from the decay of long-lived resonances and 
unstable particles (u, 77, rf , K°, A, etc.), and inaccurate 
Coulomb corrections |2|. 

A better determination of the source chaoticity is pos- 
sible by using three-particle correlations. Normalizing 
the three-pion correlation function appropriately by the 
two-pion correlator, the effects from particle misidcn- 
tification and decay contributions can be removed 
thereby isolating possible coherence effects in the parti- 
cle emission process. The resulting three-pion correlator 
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r 3 provides the means of extracting the degree of source 
chaoticity by examining its value in the limit of zero rel- 
ative momentum. Recent measurements at the CERN 
SPS from experiments NA44 and WA98 have focused on 
extracting r 3 from three-pion correlations 0, . While 
these studies have produced results which are consistent 
with a chaotic source for Pb+Pb collisions (^/snn = 17 
GeV), NA44 in particular has shown for S+Pb collisions 
(■Jsnn — 20 GeV) a result which does not appear to 
be consistent with the chaotic assumption. All of these 
prior results suffer from low statistics which limits their 
significance. We present here using charged pions the 
first high-statistics heavy-ion study of three-pion corre- 
lations, resulting in the first accurate measurement of 
the degree of chaoticity in Au+Au collisions at the Rel- 
ativistic Heavy Ion Collider (RHIC). Note that a similar 
study has recently been carried out for CERN LEP e + e~ 
collisions |(| which reports a fully chaotic source. 



J 



The present three-pion correlation study by the STAR 
experiment at RHIC supplements the published two-pion 
correlation data from Au+Au collisions at ^/snJi = 130 
GeV • A summary of the three-pion results will be pre- 
sented for two multiplicity classes. By looking at collision 
classes with different multiplicities we can vary the im- 
pact parameter, and thus the number of initially colliding 
nucleons, and study the effect of the size of the colliding 
system on the source chaoticity. We discuss the method 
of normalization of the correlation function and its ex- 
trapolation to vanishing relative momentum in order to 
extract the source chaoticity; we estimate the various sys- 
tematic uncertainties associated with these procedures. 

Before presenting our experimental results, we first 
outline the formalism which guided our analysis (for de- 
tails see Ref. 0,0 and references therein). The measured 
observable is the normalized three-pion correlator: 



r 3 (Qa) = 



(C 3 (Q 3 ) - 1) - (C 2 (Qia) - 1) - (C 2 (Q 23 ) - 1) - (C 2 (Qai) - 1) 



V(c 2 (Qi 2 ) 



1) (C 2 (Q 23 ) 

I 



1) (C a (Q31) - 1) 



(1) 



Here Q 3 = \JQ\ 2 + Ql 3 + and = J-foi-Pj) 2 
are the standard invariant relative momenta |j, |5| which 
can be computed for each pion triplet from the three mea- 
sured momenta (p l7 p 2 ,p 3 ). C 2 (p t ,Pj) = p/^pf^) = 

C 2 (Q y ) and C 3 ( Pl ,p 2 ,p 3 ) = Pl( Ztill'k( P J = 
C 3 (Q 3 ), where P represents the momentum probabil- 
ity distribution. In Ref. Q the ratio r 3 is defined in 
terms of functions which depend on all 9 components of 
(Pi>P 2 'P3); however, limited statistics even in our high- 
statistics sample requires a projection of both the numer- 
ator and denominator onto a single momentum variable, 

Qa- 

For fully chaotic sources r 3 /2 approaches unity as all 
relative momenta (and thus Q 3 ) go to zero. If the source 
is partially coherent, a relationship can be established Q 
between the limiting value of the three-pion correlator at 
Q 3 = and e, the fraction of pions which are emitted 
chaotically from the pion source (0 < e < 1): 



1 in n\ r 3 ~ 2g 
~r 3 (Q 3 =0) = ^/i^-^. 



(2) 



e gives an upper limit on the value of the two-pion A- 
parameter, which is sensitive to the fraction of coherent 
pairs in a sample, i.e. A = e(2 — e) assuming no other 
effects on A such as long-lived resonances 2J . Eq. (2) 
is not affected by the projection onto a single relative 
momentum variable. To exploit it and extract the degree 
of chaoticity e, the measured data for r 3 must, however, 
be extrapolated from finite Q 3 to Q 3 = 0. 



Similar to the two-boson correlation function, the 
three-boson correlation function C 3 (Q 3 ) is calculated 
from the data by taking the ratio gT^f] arid normalizing 

it to unity at large Q 3 . Here A (Q 3 ) = is the three- 
pion distribution as a function of the invariant three-pion 
relative momentum, integrated over the total momentum 
of the pion triplet as well as all other relative momentum 
components. It is obtained by taking three pions from a 
single event, calculating Q 3 , and binning the results in a 
histogram. B (Q 3 ) is the analogous mixed event distri- 
bution which is computed by taking a single pion from 
each of three separate events. Because of the zero in 
the denominator of the normalized three-pion correlator 
r 3 at large Qij 1 the particular method of normalization 
of C 2 and C 3 can have a strong effect on the calcula- 
tion. The propagation of statistical errors through the r 3 
functions, however, accounts for these effects completely. 
In fact, it is only with the very high statistics available 
from STAR that the calculation can be considered in the 
range 15 < Q 3 < 120MeV/c. This range is large enough 
to provide reliable extrapolation to Q 3 = 0. 

Data for the present results are from about 300K 
events taken during the ^snn = 130 GeV Au+Au run 
at STAR using the Time Projection Chamber (TPC) 
as the primary tracking detector. In the discussion that 
follows, all phase space cuts and experimental correc- 
tions are similar to the two-pion HBT analysis Q. A 
set of multiplicity classes was created by taking the 12% 
most central for the high-multiplicity class and the next 
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FIG. 1: Three-pion correlation function for central Au- 
Au events using ir~ triplets. Statistical and statisti- 
cal+systematic errors are shown. 



20% most central for the mid-multiplicity class. For both 
multiplicity bins, tracks were constrained to have px in 
the range 0.125 < pr < 0.5GeV/c, and pseudorapidity 
\rj\ < 1.0. A vertex cut was also applied to events such 
that the vertex along the z-axis (beam direction) had to 
fall within ±75 cm of the center of the detector. In the 
range 15 < Q3 < 120MeV/c, approximately 150 million 
triplets were included in both the negative and positive 
pion studies. 

The Ci correlation function was corrected for Coulomb 
repulsion with a finite Gaussian source approximation, 
using an integration of Coulomb wave functions [lOj. 
In calculating C3, the correction was applied by taking 
the product of three two-pion corrections, obtained from 
the three possible pairs formed from each mixed-event 
triplet. This type of correction approximates the three- 
body Coulomb problem to first order 0, [^. Other 
methods to more accurately estimate the true three- body 
Coulomb effect show a 5-10% smaller correction |l3j . 
This difference was applied to the Coulomb correction 
factor in calculating C3 , and the resulting shifts in the 7-3 
function were found to be within systematic uncertain- 
ties. A separate study examined the effect of inappro- 
priately applying the Coulomb correction to pions which 
come from long-lived resonances [lij . Using a rescatter- 
ing model [lflj . the value of r% was found to increase by 
10% when pairs and triplets of pions which contain pions 
from long-lived resonances were inappropriately Coulomb 
corrected. Effects of finite momentum resolution on r% 
were also studied using this model and were found to be 
insignificant. The la uncertainty in determining Q3 is 
found to be about 10 MeV/c. 

Figure shows the C3 correlation function for nega- 
tively charged pions in the high-multiplicity bin. The 
shape of C3 is mostly built up of products of two-pion 
correlations with the effect of true three-pion correlations 
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FIG. 2: Calculation of rz for (a) central and (b) mid-central 
ir~ events, and (c) central and (d) mid-central 7r + events. The 
fits shown use Eq. Statistical and statistical+systematic 
uncertainties are shown. 



being more subtle. At large Q3> C3 approaches unity and 
for an ideal pion source, i.e. A = 1, C3 would approach 
6 at Q3 = (this is not the present case since A < 1). 
A Gaussian parameterization is inadequate to describe 
this correlation function; this is consistent with results 
obtained in other experiments and a 

simulation[lHE3. 

In calculating 7*3, the actual binned values of the cor- 
relation function for the various values of Q3 are used 
instead of a fit 15]. In order to use Eq.([TJ, triplets are 
obtained that pass all of the momentum space and ex- 
perimental cuts. Q3, Q12, Q23 and Q31 are calculated 
from the triplet and the three pairs that can be formed 
from the triplet. The values C 3 (Q3), C 2 (Q12), C% (Q23) 
and C2 (Q31) are then computed from the binned actual 
two- and three-pion correlation functions. These values 
are then used to calculate 7-3, which is then binned as a 
function of Q3. The average for each bin is then calcu- 
lated to obtain the final result. Systematic uncertainties 
are greatest at the low Q3 end due to track merging ef- 
fects and the uncertainty in the Coulomb correction. The 
parameters controlling these effects were modified ±20% 
from the nominal values to obtain the overall systematic 
uncertainty in each bin. 

The results for the two multiplicity bins are shown in 
Figure[2]for 7r~ and 7r + , plotted as functions of Q\. Plot- 
ting in this way is suggested by the theoretical analysis 
in [3j which shows that the leading relative momentum 
dependencies in the numerator and denominator of Eq. 
(|T|l are quadratic [Hi ] , allowing for a linear extrapolation 
of the results shown in Figure[2]to Q3 — by fitting them 
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to the form 

r 3 (Q 3 )/2 = r 3 (0)/2-aQl (3) 

where r 3 (0)/2 and a are fit parameters. From Figure 
121 it appears that the normalized three-pion correlator 
r 3 (Q 3 ) does indeed show a leading quadratic dependence 
for the smaller Q§ values (Eq. © was fit to the range 
< Q 3 < 60 MeV/c). 

The resulting intercepts r 3 (0)/2 are shown in Figure |3J 
along with the results of WA98 and NA44. Error bars for 
STAR points are statistical+systematic. As mentioned 
earlier, the systematic error is computed by varying sev- 
eral parameters independently including particle track 
cut parameters. The variation of the parameters is seen 
to produce, in general, asymmetric variations in the ex- 
tracted intercepts. Intercepts from the quadratic fits as 
well as from quartic fits (i.e. adding a quartic term to 
Eq. (3) and fitting over the broader range < Q 3 < 120 
MeV/c) are shown for comparison, and are seen to agree 
within errors. The STAR ir + and ir~ results are also 
seen to agree within error bars. NA44 reported a result 
close to unity for Pb-Pb interactions, but a much lower 
result for S-Pb Q, both with no clear Q 3 dependence. 
The Pb-Pb result from WA98 is somewhat smaller than 
that from NA44, although they agree within error bars, 
and the ^-dependence in their result is similar to what 
is seen in STAR Q. 

Figure |U shows the results from the calculation of e for 
STAR's measurements, and for those from WA98 and 
NA44, plotted versus charged particle multiplicity per 
unit pseudorapidity, dN/dr). The calculation was done 
starting with the results of Figure |31 decreasing them 
by 10% to approximately take into account the overcor- 
rection produced by Coulomb-correcting long-lived reso- 
nances (see earlier discussion) and using Eq. (2). It was 
assumed that the 10% correction also applies to the SPS 
data, and to be conservative, a ±5% systematic uncer- 
tainty on the correction (i.e. 10% ± 5%) was included in 
all of the error bars shown. The plot shows an increasing 
trend in the STAR 7r~ and 7r + results going from mid- 
central to central collisions. For the mid-central data, the 
results for e show a partially chaotic source, as seen in 
the SPS results. The central data appear to give a mostly 
chaotic pion source. Including the SPS measurements 
into the overall systematics, there appears to be, within 
the uncertainties shown, a systematic increase in s with 
increasing dN/drj, the smallest value being for SPS S-Pb 
collisions and the largest value for STAR central Au-Au 
collisions {dN/drj for charged particles at mid-77 for SPS 
S-Pb, SPS Pb-Pb, STAR mid-central, and STAR central 
are approximately 100(scaled from S-S), 370, 280, and 
510, respectively [H El)- It is also found for the STAR 
results that the upper limit on the two-pion A-parametcr 
obtained from e using the relationship mentioned earlier 
is in the range 0.71 — 0.81 for mid-central and 0.91 — 0.97 
for central events. The actual values for A extracted from 
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FIG. 3: rs(Q3 = 0)/2 from STAR and two other experiments 
0,0. For STAR, (a) central and (b) mid-central results are 
shown for 7r~ (circles) and ir + (squares) data. The other 
experiments use n~ data only. STAR results for fitting with 
both a quadratic and quartic functions are shown. 

STAR it- - ir~ HBT measurements are 0.53±0.02 for 
mid-central and 0.50±0.01 for central events (the ir + —ir + 
values agree with these within errors) 7J. The lower A- 
values extracted from the two-pion experiment can be 
explained in terms of long-lived resonance effects, which 
nicely cancel out in a three-pion analysis |l5|. 

In summary, we have presented three-pion HBT results 
for y/SNN = 130 GeV data at STAR, and have shown that 
for the central multiplicity class the STAR data indi- 
cate a large degree of chaoticity in the source at freeze- 
out, whereas for the mid-central class the source is less 
chaotic. Our r 3 results are close to those extracted in 
SPS Pb+Pb collisions, but differ from the low value ob- 
tained in SPS S+Pb collisions. The comparison between 
SPS and STAR results suggests a systematic dependence 
of the chaoticity on particle multiplicity. High statistics 
from STAR have allowed a normalized three-pion corre- 
lator calculation that extends to 120 MeV/c in Q 3 , and 
the dependence on this variable has been shown to be 
quadratic in nature for low Q 3 . STAR's measured val- 
ues provide increased confidence in the validity of stan- 
dard HBT analyses based on the assumption of a chaotic 
source for central collisions at RHIC. 
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